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Abstract. We consider the problem of resolving all pairwise interactions of shock waves, con- 
tact waves, and rarefaction waves in 1-dimensional flow of an ideal polytropic gas. Resolving 
an interaction means here to determine the types of the three outgoing (backward, contact, and 
forward) waves in the Riemann problem defined by the extreme left and right states of the two 
incoming waves, together with possible vacuum formation. This problem has been considered by 
several authors and turns out to be surprisingly involved. For each type of interaction (head-on, 
involving a contact, or overtaking) the outcome depends on the strengths of the incoming waves. 
In the case of overtaking waves the type of the reflected wave also depends on the value of the 
adiabatic constant. Our analysis provides a complete breakdown and gives the exact outcome of 
each interaction. 
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1. Introduction 

We consider the problem of resolving all interactions of pairs of elementary waves in one-dimensional 



compressible, adiabatic flow of an ideal and polytropic gas as described by the Euler system (2.1 1, 



(2.2), (2.3). The elementary waves are shocks, contacts, and centered rarefactions. The interactions 
are treated as in the Glimm scheme [3j: if the incoming waves connect a left state £/;, via a middle 
state U m , to right state U r , then the "interaction problem" is to determine the type and strength of 
the resulting outgoing waves obtained by resolving the Riemann problem (Ui,U r ). Here U denotes 
any triple needed to specify the state the gas; we will mostly work with specific volume r = l/p 
(p = density), fluid velocity u, and pressure p. 

For interactions that involve only shocks and/or contacts the solution of the Riemann problem 
(17;, U r ) provides the exact solution of the wave interaction. When one of the waves is a rarefaction 
the actual interaction involves "penetration" of a rarefaction wave by a shock, a contact or by 
another rarefaction. In these cases the exact solution is more involved. We do not treat the problem 
of penetration and consider only the Riemann problem defined by the extreme states Ui and U r . 

Remark 1.1. In cases of wave penetration the extreme Riemann problem (Ui, U r ) may not provide 
an accurate description of the asymptotic behavior of the wave interaction. E.g., in an overtaking 
interaction of a shock and a rarefaction, the shock may not pass through the rarefaction ( incomplete 
penetration), and in this case the solution of the Riemann problem (Ui,U r ) presumably does not 
describe the asymptotic behavior in the exact solution. See Part D of Chapter III in j2], Section 
3.5.2 in [I], and [4| for further details and references. 

Many authors have studied the problem of pairwise interactions. Since it gives rise to Riemann 
problems it may be said to originate with Riemann's and Hugoniot's pioneering works [5j|6|ll] on gas 
dynamics. According to [2], Jouguet [7] in his work on combustion considered interaction phenomena. 
A systematic approach was initiated by von Neumann [9] in his work on hydrodynamics. Among 
other issues he considered the interaction of two shock waves, and called attention to the fact that a 
contact wave emerges in such interactions. He also observed that in interactions of overtaking shocks 
in an ideal gas, the reflected wave is necessarily a rarefaction whenever the adiabatic constant 7 
satisfies 7 < f . Shock-shock and shock-contact interactions were analyzed in detail by Rozdestvenskh 



& Janenk o [12] ; see also Smoller 10 . A comprehensive treatment is presented by Courant and 
Friedrichs |2] who also list results for shock-contact and rarefaction-contact interactions. A partial 
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treatment is also given in |8| |12| . The most complete work on pairwise interactions in ideal polytropic 
gases to date appears to be the analysis by Chang and Hsiao [l]. These authors considered all the 



ten essentially different types of pairwise interactions (see (3.1 1-(3.3|), for all values of the adiabatic 
constant 7 > 1. In particular they list the possible outcomes of all ten interactions. However their 
analysis does not completely delimit all the various outcomes. 

In revisiting the interaction problem our objectives are to provide a complete breakdown of all 
possible cases, and to determine the exact outcome of each interaction. In addition we specify 
precisely when a vacuum state appears among the outgoing waves. Our analysis is motivated by a 
desire to investigate relevant measures for wave-strengths in large data solutions of the Euler system. 
The latter issue will be pursued elsewhere. 

Let us explain what we mean by "complete breakdown of all cases." This requires a little back- 
ground. We work in a Lagrangian frame such that contact waves appear stationary, and the primary 
unknowns are specific volume r, particle velocity u, and pressure p (alternatively, specific entropy 
S). We parametrize shocks and rarefactions by pressure ratios across the wave, while contacts are 
parametrized by specific volume ratios. With a slight abuse of terminology "a shock /," say, refers 
to a shock wave connecting two states U- (to the left of the shock) and U+ (to the right of the 
shock), and such that the pressures satisfy ^j^y = /• We also refer to / as the strength of the 
wave. A wave changes type (shock <-> rarefaction, or up-contact ■<-> down-contact) as this pressure 
(or specific volume) ratio crosses the value 1. 

Now, an interaction of two incoming waves x and y generically gives rise to three outgoing waves: 
a backward wave B, a contact C, and a forward wave F, and, possibly, a vacuum state. The 
interaction problem amounts to determining whether a vacuum occurs, together with the values of 
B, C, F, in terms of x, y. 

It does not seem possible to determine the functions B(x,y), C(x,y), and F(x,y) explicitly 
in all cases. A somewhat weaker request is to ask for the "transitional curves" in the plane of 
incoming strengths. That is, to determine the curves in the (x, y)-plane which delimit the regions 
{(x, y) I B(x, y) ^ 1}, {(x, y) \ C(x, y) ^ 1}, {(x, y) \ F(x, y) ^ 1}, as well as the vacuums regions. (It 
turns out that these transitional sets are either C 2 -curves, or finite unions of such.) It is possible to 
find explicit equations for these curves and to determine their relative locations and intersections, 
and thereby determine all possible outcomes in all pairwise interactions. 

In this work we carry out the detailed computations that yield a complete breakdown in this 
sense. In doing so we improve on the existing results in the literature. E.g., one of the more involved 
cases occurs when a backward shock is overtaken by a backward rarefaction (S R) in a gas with 
adiabatic exponent 7 < |. Depending on the strengths x, y of the incoming waves, the interaction 
may produce any one of four outgoing wave configurations: 

SJS, RJS, SJR, or RJR, 

with possible vacuum formation in the latter case, see Figure 5 (left diagram). Here R, J, and S 
denote rarefactions, contacts, and shocks, respectively; an arrow indicates the direction the wave 
is moving relative to the Lagrangian frame, and ^ indicate increase/decrease of density across a 
contact. In our analysis we delimit the various cases in terms of certain curves in the (x, y)-plane. 
These curves are given in terms of explicit equations - in some cases as explicit graphs y(x) or x{y). 

The calculations are surprisingly involved, with a number of sub-cases for various combinations 
of strengths x, y and values of 7. It is possible that there are more efficient ways of parametrizing 
waves than the present choice of pressure and specific volume ratios. In fact, for some cases we find 
it useful to consider the entropy variable as well. However, we have not been able to determine more 
efficient parameters that would, e.g., provide explicit formulas for outgoing strengths in terms of 
incoming strengths. 

The rest of the article is organized as follows. In Section [2] we recall the Euler system and two 
parametrizations of its wave curves, one in (t, u,p)-space and the other in (r, u, S')-space. In Section 
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|2.3| we briefly review the solution of the Riemann problem and characterize the occurrence of a 
vacuum. We also list and group the ten essentially different pairwise interactions, modulo reflections 
of the spatial variable x O —x. We analyze head-on (Group I) interactions in Section |4j and 
interactions with a contact (Group II) in Section [5[ The results for these groups are summed up 
in Theorems |3.2| and |3.3[ All of the remaining analysis deals with the more involved cases where a 
shock or rarefaction overtakes another shock or rarefaction (Group III). The results for overtaking 
interactions are given in Theorem |3.4| It turns out that the outcome depends not only on the 
incoming wave-strengths but also on the adiabatic constant 7. Considering (for concreteness) the 
case of overtaking backward waves we divide the analysis into three sub-groups: S S (Ilia), SR 
(111b), and RS (Hie) interactions. As observed by Chang and Hsiao Ej there are two values (7 = | 
and 7 = 2) at which the global structure of the transitional curve for the reflected wave (i.e., the 
curve {(x,y) | F(x,y) — I}) changes character. The type of the transmitted (backward) wave for 
Group III interactions is analyzed in Section [6j In Section [7] we treat the reflected wave; this is 
where the most involved analysis occurs, due to the different behaviors for 7 ^ | and 7^2. To 
analyze the outgoing contact it turns out to be convenient to use (r, u, S) variables; the details are 
given in Section [8] Finally, in Section [9] we analyze vacuum formation in overtaking interactions. 
Section [TO] collects various auxiliary results needed in the analysis. 

2. The Euler system, wave curves, Riemann problems, and vacuum 

2.1. The 1-d compressible Euler system for an ideal gas. We consider 1-d compressible flow 
of an ideal, polytropic gas. The conservation laws for mass, momentum, and energy are given by 

T t -u x = (2.1) 
u t +p x = (2.2) 
E t + (up) x = 0, (2.3) 

where x denoted 1 a Lagrangian (mass) coordinate, r = 1/p is specific volume, u is the fluid velocity, 

2 

p is pressure, and E = e + ^ is specific total energy. Here e denotes specific internal energy, which 
for an ideal gas is given by 

pr 

e(r, p) = , where 7 > 1 is the adiabatic constant. (2-4) 

7-I 

For later reference we note that the local sound speed c is given by 

c 2 =7rp. (2.5) 

2.2. Wave curves. Each fixed state (r, U,p) G M. + x M x M+ has three associated wave curves 
consisting of those states that can be connected to (f, u,p) (on its right) by a single backward wave, 
a single contact, or a single forward wave, respectively. The backward and forward waves arc cither 
rarefactions or entropy admissible (i.e. compressive) shocks. We use pressure ratios Pright/Pleft to 
parametrize the wave curves in the extremal fields, and specific volume ratios T r i g ht/Ti e ft for the 
contact field. We let b, c, / denote the pressure ratio (strength) of a backward, contact, and forward 
wave, respectively. Using the expressions from [10] (p. 354) the wave curves are given as follows: 

WUr.r. „.,,)= I u-r(l,h'TT> I I 5 ? ? < f < 1 (2-6) 



S : b>l 



J(c; f,u,p)=[ u\ l^c (c>0) (2.7) 




and t will later be used for other parameters as well. The context will make it clear which use is intended. 



PAIRWISE WAVE INTERACTIONS IN IDEAL POLYTROPIC GASES 



5 



<f>(f)f 

W{f;r,u,p) = ( u+i>{f)jf§ 
fP 

The auxiliary functions <j) 1 <j>, ip, and ip are given by 



B 

S 



/>1 

o</< 1 



(2.8) 



0(6) 



'(/) 



/ 



1+q6 
6+a 

-1/7 



1+a/ 




- 1) 



«(&-!) 
s/b+a 



«(/-!) 



0</<l 

where the constants ft, i/, £ are given in terms of the parameter^] 

7-1 




7 + 1 



as follows: 



K = Vi-a e (0,1) , 



vT 



€ (0,1) 



e (O.+oo), C 



(2.9) 



l + o 



6(0 A). 



We note that a is a strictly increasing function of 7, and that 



^> 7 



and 



7 



10 



The auxiliary functions cf>, (j) , ?/> , and ?/> are analyzed in Section 

It will be convenient to perform some of the calculations in (r, u, £) -variables (see Section |8|, 
where S — specific entropy, and we proceed to record the expressions for the wave curves in these 
variables. The pressure is now given by 



p = p(t, S) = Kt 7 exp (jf^j {K, c v positive constants). 



(2.10) 



We parametrize all wave curves in (r, u, S')-space by specific volume ratios Trightl^l&ft- We use I for 
backward waves, c for contacts (as above), and 1 for forward waves. From the parametrizations in 
(t, it, p)-space, together with (2.10), we obtain the following expressions for the wave curves (r, u, S)- 
variables: 

If 

§ + c v rf(l) 




R 

S 



I > 1 
a < I < 1 




1 § c (c> 0) 



W(i; t, u, S) = u+^(b)^/fp 
\ S + c v rj(i.) 

The auxiliary function £ , £ , 77 , rj are given by 



R 
S 



< L < 1 
1< i < i 



(2.11) 
(2.12) 
(2.13) 



£(0 



/>i 



Vo+ 1 



(i-O 

v i— a 



77(0 



< 1 < 1 



, /P(l-aZ)\ 



i > 1 

a < £ < 1 



^For brevity we use a instead of the more common notation fi 2 . 
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y/a+1 



(1-0 




7?(t) 



< L < 1 



1 < K 



Remark 2.1. j4s is clear from the parametrizations (2.11)-(2.13) the entropy remains constant 
across rarefactions, while it necessarily jumps across contacts. 



2.3. Riemann problems and vacuum formation. A Riemann problem refers to the problem of 
resolving a single, initial jump discontinuity connecting two constant states. It is known that the 
Euler system for an ideal, polytropic gas admits a self-similar solution to any Riemann problem. 
Futhermore, this solution is unique within the class of self-similar solutions consisting of constant 
states connected by compressible shocks, contacts, and centered rarefactions, possibly including a 
vacuum state; see (2j. 



Consider the Riemann problem for (2.1)-(2.3) with left state (f,u,p) and right state (r,u,p). We 
use capital letters B, C, and F to denote the resulting outgoing backward, contact and forward 
waves, respectively. Traversing the resulting wave- fan from left to right and using the expressions 



fl2.6M2.8D, yield 



= ip(F)yC B 4>(B) — ip(B) , p = BFp. 



(2.14) 



By using the first and last equations in the middle equation, together with the relation (10.1), we 
obtain the following equation for B: 

'pB" 
P 



F(B;r,p,f,p)=F(B). 



(2.15) 



\/Tp y Tp 

To solve the Riemann problem ((f,u,p), (r, u,p)\ amounts to determining the pressure ratio B of 
the outgoing backward wave from J-(B) = -^SS. ^From this one then determines C and F from 



(2.14) 



The function ip is strictly increasing and tends to +00 at +00 (see Section 10). It follows that 
the map B t— > J-(B; T,p, f,p) has the same properties, and that the Riemann problem has a unique 
solution without vacuum if and only if (f, u,p) and (r, u,p) are such that 

-Tin - -\ u — u 

■F(0;t,p,t,p) < 



'Tp 



This is the case if and only if 



u < 



7-1 



(c + c) 



(no vaccum), 



(2.16) 



where the sound speeds c and c are given by (2.5). If the condition (2.16) is not satisfied, then 



we agree to solve the Riemann problem as follows: a backward rarefaction connects (f, u,p) to 
(oo,0,0), followed by a forward rarefaction connecting (00, 0,0) to (t,u,p). In the (x,i)-plane the 
two rarefactions are separated by a vertical line along which r = 00. The two rarefactions, together 
with the line x = 0, span the fan (see 10 




Remark 2.2. We note that the rarefactions on either side of a vacuum may be in different entropy 
states. In particular this can occur in vacuums resulting from interactions: see region {c < 1 < /} 



in Figure 3, and region {y < 1 < x} in Figures 4-6. Recalling Remark 2.1 we denote the outcome 



of such interactions by RJR according to whether the left and right entropy states satisfy Si S r . 
Thus, for these outcomes, J indicates the direction of the entropy jump across the vacuum. 
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Figure 1. Interactions: head-on (I), involving a contact (II), overtaking (III). (Schematic) 



3. PAIRWISE INTERACTIONS AND SUMMARY OF RESULTS 

3.1. Pairwise interactions. We shall consider all possible wave configurations where two elemen- 
tary waves (compressible shocks, contacts, and centered rarefactions) interact. We collect the inter- 
actions into three groups: head-on (I), involving a contact (II), and overtaking (III), see Figure 1. In 
interactions, lower case letters are used for the incoming waves (/ for forward waves, c for contacts, 
and b for backward waves), while capital letters B, C, F, denote outgoing waves. In Figure 1 the 
waves are depicted schematically as if no penetration occurs: all waves, including possible incoming 
and outgoing rarefaction fans, are drawn as single lines in the (x, t)-plane. As described in the 
introduction the interaction problem is to resolve the particular Riemann problem defined by the 
extreme left and right states in the incoming waves. 

Due to the symmetry x o —x there are ten essentially different pairwise interactions: 



la : 


SS 


lb : 


SR 


Ic : 


RR 




_a > 


Ila : 


S J 




_a < 


lib : 


SJ 




_k > 


lie : 


RJ 




_^ < 


lid : 


RJ 


Ilia : 


SS 


Illb : 


SR 


IIIc : 


RS 



head-on interactions 



(3.1) 



interactions with a contact 



overtaking interactions. 



(3.2) 



(3.3) 



Each interaction defines a particular Riemann problem, the interaction Riemann problem, for which 
we obtain a non-linear algebraic equation determining the strength of the outgoing backward wave, 
say, as detailed above for a general Riemann problem. When expressed in terms of the auxiliary 
functions 4> , <fi, tp , ip, this nonlinear equation is the same within each of the groups I - III. However, 
a breakdown into individual cases is necessary to determine exactly the outcomes of each type of 
interaction. 



Remark 3.1. The parametrizations of the wave curves in Section 2.2 are in terms of pressure and 
specific volume ratios. This has the following useful consequence: the outgoing wave strengths in 
interactions are given in terms of three nonlinear algebraic equations that involves only the outgoing 
wave strengths (and j). In particular, the equations are independent of say, the left-most state in 
the interaction. See equations (4.1)-(4.3), (5.1)-(5.3), and (6.1)-(6.3). 



3.2. Summary of results. In this section we collect our results for interactions in each of the three 
groups. For each group we consider every combination of incoming waves, as listed in (3.1)-(3.3), and 
we describe all possible combinations of outgoing waves in the corresponding interaction Riemann 
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problem. As discussed in the introduction it does not seem possible to determine outgoing strengths 
as explicit functions of incoming strengths. Instead we describe the "transitional curves" in the plane 
of incoming strengths. For example, for head-on interactions with incoming wave strengths b and /, 



the outgoing backward wave B = B(b, f) may be a rarefaction or a shock. According to (2.6) these 
outcomes correspond to B(b,f) < 1 and B(b, f) > 1, respectively. The transitional curve for the 
backward wave in this case is the curve {(b, f) | B(b, /) = 1}. Similarly, there are transitional curves 
for the outgoing contact C and forward wave F. It will turn out that these sets really are curves in 
the plane of incoming strengths. In addition we determine the location of all "vacuum-curves", i.e., 
the curves in the plane of incoming strengths which delimit interactions for which the interaction 
Riemann problem contains a vacuum state. 

To describe the results we make reference to Figures 2-6. First, for Group I and Group II 
we give two diagrams in the plane of incoming strengths (i.e., the (6, /)-plane for Group I and the 
(/, c)-plane for Group II). The left diagram depicts the various combinations of incoming waves: see 
left diagrams in Figure 2 and Figure 3. The diagram on the right in each of these figures shows the 
resulting, outgoing wave configurations, together with the transitional curves and vacuum curves. 
It turns out that the vacuum curve in both of the right diagrams in Figure 2 and Figure 3 are 
explicitly given as graphs. The same is true for the transitional curve for the outgoing contact in la 
(S S) interactions (it is given by bf = 1). However, the transitional curve for the outgoing contact 
< 

in lid (RJ) interactions is only implicitly determined in our analysis. The figures are schematic; 
see Remark 15.21 

Before considering Group III interactions we state the results for Groups I and II. We recall that 
the strength of a shock and rarefaction is defined as the pressure ratio Prightjpieft-, while the strength 
of a contact wave is defined as the specific volume ratio T r i g ht I Tie ft = Pie ft / Pright ■ As a consequence 
all incoming strengths in pairwise interactions are assumed to be ^ 1. 

Theorem 3.2 (Group I interactions). Consider the head-on interactions of two elementary waves 



in (3.1). Let the incoming forward and backward waves have strengths f and b, respectively (see left 
diagram in Figure 1). Then: 

(i) the outgoing backward wave B satisfies 

B^l & 6^1 (3.4) 

(ii) the outgoing forward wave F satisfies 

F%1 «> f%l (3.5) 

(iii) the outgoing contact C satisfies 

- Ia (SS): C^l & bf^l 

- Ib (SR): C > 1 

- Ic (RR): C = 1 (no outgoing contact) 

(iv) for vacuum formation we have: there is never vacuum formation in Ia and Ib interactions, 
while for Ic interactions vacuum occurs if and only if the incoming strengths f and b satisfy 

b c + ,r c <l where C=^- (3-6) 

The situation is summarized in Figure 2. 

The proof of Theorem 3.2 is given in Section [4] Some additional information about the outgoing 



strengths in Ia and Ib interactions are recorded in Lemma 4.1 

Theorem 3.3 (Group II interactions). Consider the interactions of an elementary forward wave 
with a contact wave in p?j| ). Let the incoming forward and contact waves have strengths f and c, 
respectively (see middle diagram in Figure 1). Then: 
(i) the outgoing backward wave B satisfies 

- Ha and lib (SJ): B^l <=> c ^ 1 . (3.7) 
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Ic: RR 



(RS) 



lb: SR 



la: SS 



vacuum 






RR / 


(SJR) 




RR 






RJS 


SJS 






SJS 


bf=l 



Figure 2. Group I interactions: incoming (left) and outgoing (right) waves; 
schematic. For vacuum configurations, see Remark |2.2| 



- lie and lid (RJ): 
(ii) the outgoing forward wave F satisfies 



B ^ 1 



F ^ 1 



/£1 



(3.8) 



(3.9) 



(iii) the outgoing contact C satisfies 

- Ha (SJ). 

- lib (SJ). 

- lie (RJ): 

- lid (RJ) 

< 



c < C < 1 
1< C < c 
C = c< 1 

in this case /, c > 1 and C < c. However, the outgoing contact may be 

> 

either a J (C > 1) or a J (C < 1) contact. More precisely, for fixed c > 1 the map 

f^C = C(f,c) 

takes on every value in (0,c) as f increases from 1 to oo. 
(iv) for vacuum formation we have: there is never vacuum formation in Ha, lib, and lid inter- 
actions, while for He interactions a vacuum occurs if and only if the incoming strengths f 
and c satisfy 

/>/*(c):=( I JL=)* where (=^. (3.10) 

The situation is summarized in Figure 3. 

The proof of Theorem 3.3 is given in Section [5] The location of the transition curve {C — 1} in 
lid interactions is analyzed in Section [5. 6| 

We finally consider Group III interactions of two overtaking waves, which we take to be backward 
waves. The possible combinations of incoming waves are S S (Ilia), SR (IHb), and RS (IIIc), since 
two backward rarefactions do not meet. Letting x and y denote the strengths of the left and right 
incoming waves, respectively, the incoming configurations are listed in the (x,y)-plane in Figure 4. 
As described in the introduction the outcomes of these interactions depend sensitively on the value 
of the adiabatic constant 7 (as well as on the incoming strengths), and we consider four different 
cases: 7 G (1, |), 7 = |, 7 £ (§>2), and 7 > 2. The resulting outgoing configurations are given in 
the (x, y)-plane by Figure 5 for 7 £ (1, |], and in Figure 6 for 7 > |. 
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lib: SI 


lid: RJ 


Ha: SI 


Hc:RJ 











SJR 










RJS 




SJR 


SJS 


RJR 


RJR 






vacuum 



1 f*(0) 



Figure 3. Group II interactions: incoming (left) and outgoing (right) waves; 
schematic (see Remark 5.2). For vacuum configurations, see Remark 2.2 



In the following statement y denotes the unique root different from 1 of the function a(y) 



4>(y) — il>(y) when 7 7^ f , while yo — 1 when 7 = |. (The function a is analyzed in Section 10.3 ) 



Also, to shorten notation we write e.g. B instead of B(x, y) for the strength of the backward outgoing 
wave. 

Theorem 3.4 (Group III interactions). Consider the interactions of two overtaking, elementary 
backward waves listed in Let the left and right incoming waves have strengths x and y, 

respectively, and let the outgoing waves have strengths B, C and F (see right diagram in Figure 1). 
Then for all values 0/7 > 1 : 

(i) the locus {B = 1}, indicated by a dashed curve in Figures 5-6, coincides with a C 2 -smooth 
and decreasing graph y — k(x). The graph satisfies k(l) — 1 and has asymptotes x = and 



y — y{l) £ (0,1) where 2/(7) is given in (9.5). The outgoing backward wave B is a shock 
(i.e. B > 1) if and only if y > k{x) . 
(ii) the locus {F = 1} consists of three curves, indicated by solid lines in Figures 5 - 6: {x = 1}, 
{y = 1}, and {y — h(x)}. We refer to these figures for the type of the outgoing forward 
(reflected) wave in Group III interactions (see also Figures 8-14). As indicated there the 
graph y — h(x) have different global properties according to the value of the adiabatic constant 
7 > 1. However, for each value 0/7 > 1; 

— y = h(x) has a horizontal asymptote y* = y«(~f) as x f +oo, 

— the graphs y = k{x) and y — h(x) intersect at the point (i, y n ). 
(hi) the outgoing contact C satisfies: 

— C < 1 in S S interactions 

— C > 1 in S R and RS interactions. 

(iv) for vacuum formation we have: a vacuum never appears in S S and RS interactions, while 
a vacuum appears in an S R interaction if and only if the incoming strengths y < 1 < x 
satisfy 



0<y< 

where 



III— ^1 



(3.11) 

v^Jx + a + k(x — 1) 

The situation is summarized in Figures 4, 5, and 6. 



( = and v(x) .— . 

7 v x + ax 2 



PAIRWISE 
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Figure 4. Group III interactions: incoming waves. 
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no interaction 
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" vacuum RJ R 







\ sis 






SJR 


RJS \ 






XT~~~ SJS 




SJR 


no interaction 










vacuum RJ R 



l<y<5/3 



Y=5/3 



Figure 5. Group III interactions: outgoing waves when 1 < 7 < | (left) and when 



7=| (right); schematic. For vacuum configurations, see Remark 2.2 



Theorem ^4 follows from the propositions in Sections [6j[7j[8j and|9] The analysis in these sections 
provide further details about the location and intersection properties of the graphs y = k(x) and 
y = h(x); see Figures 8-14. When possible we also give explicit expressions for k(x) and h(x). 



Remark 3.5. Before starting on the proofs of Theorems ] 3. S\ \ 373[ and \3.4\ we note that the analysis 
will make use of a number of auxiliary functions and results about their properties. These are for 
the most part collected in Section [75| 



4. Proof of Theorem 13.2 



4.1. Group I: head-on interactions. Consider the interactions in Group I listed in (3.1). Re- 



ferring to the left diagram in Figure 1 we use ( 2.6 )-( 2.8 1 to traverse the waves before and after 



interaction. We obtain three equations for the outgoing strengths B, C, F in terms of the incoming 
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RJS \ 




SJR 




\SJR 
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SJS 


no interaction 


RJR 


vacuum 


RJR 



- < -\ SJS 
RJS . 






\ SJR 




SJR 


RJR \ 


SJS 








SJR 


no interaction 


RJR 


RJR 

vacuum 



5/3<y<2 



Y*2 



Figure 6. Group III interactions: outgoing waves when I < 7 < 2 (left) and when 



7 > 2 (right); schematic. For vacuum configurations, see Remark 2.2 
strengths b and /: 

0(6)0(/) - C%(B)${F) 

i>U) -4>{b)^ 'fHf) = -^{B) + i>{F)^CB%{B) 
bf = BF. 



(4.1) 

(4.2) 
(4.3) 



functions M and N, we obtain the following nonlinear equation for B = B(b, /): 

Q(B;b,f)=0, 



Eliminating F and C, using the relation (10.1) and the definitions (10.2), (10.3) of the auxiliary 

(4.4) 
(4.5) 



where 

Q(B; b, f) := i>(B) + %(§)M(b)N(f) + $ (/) - N(f)^(b) . 

We observe that since ip is strictly increasing, the map B h-> Q(B;b, /) has the same property. 
In what follows we first analyze the strengths of the outgoing backward and forward waves B and 



F and establish parts (i) and (ii) of Theorem 3.2 This information is then used to prove the claims 
about the outgoing contact C. Finally we establish the claims about vacuum formation in head-on 
interactions. 



4.2. Proof of Theorem 3.2 part (i). By the defining equation (4.4) for B — B(b,f) and the 
monotonicity of B M> Q(B;b,f), we have by (10.2) that 

B^l & g(B;bJ)=0^g(l;b,f)=N(f)g(b,f), (4.6) 



where N(f) > and 



g(b,f) := te)M(6) - ^(6) - $()) 



Thus, to verify (3.4) we need to show that 

5 (6,/) ^0 6^1. 
Since g(l, /) = it suffices to show that dbg(b, /) < for all b, f > 0. We have 

d b g{b, f) = M'(J#(i) - J^M(ft$'(£) - . 



(4.7) 
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Using the properties of M and ip we get 

d b g(b,f)<M'(b)^(^)<0 
On the other hand, if bf < 1 then 

1 



d b g(b,f)<M'(b)^(^ i 



whenever bf > 1 
M{b) 



<0, 



since m(b) < \ for all b > 0, while £(q) > h for all q > 1 (see Section 10 1. This proves (3.4) 



4.3. Proof of Theorem 

we have 



3.2 



bf 



part (ii). Using (4.3) and the fact that B H > Q(B; b, f) is increasing 

G(B;bJ)^G(bJ). 



F= -^^1 B^bf 

where 

G(b, f) := G(bf; 6, f) = ij(bf) + i>(f) - N(f)^(b) . 
We have G(b, 1) = 0, while 

d b G(b, f) = M'(bf) - N(f)ij'(b) , d f G(b, f) = bi>'{bf) + 4>'{f) - N'{f)^(b) . 
Using the properties of the auxiliary functions we see that each term in dfG(b, f) is non-negative 



whenever < b < 1, and that their sum in this case is strictly positive. Thus, (3.5) holds for 



< b < 1. For b > 1 we consider instead how G changes along hyperbolas bf = const. Differentiating 
in the direction of increasing / we have, for b > 1, that 



(-b,f)-V {bJ) G(b,f) = b^'(b)N(f) + M\f)- fN'(fW(b) 
> b^'(b)N(f) - fN'(f)i)(b) 
= ^(b)N(f)[£(b)-n(f)}>0 



(4.8) 



by the properties of the auxiliary functions £ and n recorded in Section 10 Again, since G(b, 1) = 0, 
it follows that (3.5) holds in the set (b, f) € {6 > 1, bf > 1}. We finally observe that 



G(l,f) = V (/) + </>(/)< whenever f<l 



which, together with (4.8), verifies (|3.5[) also in {b >l,bf< 1}. This proves (3.5). 



4.4. Proof of Theorem 3.2 part (iii). Concerning the outgoing con tact C in Group I interactions 
we first observe that C = 1 in RR (Ic) interactions. Indeed, by Remark 2.1 if C ^ 1 then the entropy 



would take different values on either side of the outgoing contact. However, we have already shown 
that the outgoing backward and forward waves are both rarefactions in this case, whence (Remark 
2.1) the entropy takes the same value on each side of the outgoing contact. 



For the other two interactions in Group I we shall use (4.1 ) together with the following 



Lemma 4.1. For la and lb interactions with incoming strengths b, f , and outgoing strengths B , C , 
F ( see Figure 1 ), the following holds: 

• For la (S S) intersections: f < F < 1 < B < b and B + F < b + f . 

• For lb (SR) interactions: F < f < 1. 

Proof. Consider first la interactions for which / < 1 < b. The properties of the auxiliary functions 
show that 

G(b; b, f) = ^(b) [1 - N(f)] + VK/) [1 - M(b)} > - g(B; b, f). 

As £/(•; b, f) is increasing this shows that b > B, and since bf — BF it follows that F > /. Combining 
this with (pSl} and (pTa), we obtain that f < F <l < B <b. But then b-f>B-F>0, whence 



(b - f) 2 > (B - F) 2 . Since bf = BF it then follows that (b + f) 2 > (B + F) 2 , i.e. B + F <b + f. 



14 



CENG CHEN, ERIK E. ENDRES, AND HELGE KRISTIAN JENSSEN 



Next, for lb interactions the incoming strengths satisfy b, f < 1. As bf — BF, we have / > F if 
and only if B > b, which is the case if and only if 

Q(B; bj) = 0> Q(b; b, f) = ^(6) [l - N(f)] + 4>{}) [l - M{b)] . 

For b, f < 1 both terms in the expression on the right are negative. This shows that / > F in lb 
interactions. □ 



to rewrite (4.1 ) as 



For a la interaction we use / < 1 < b in the explicit expressions for 4> and <f>, together with (|4.3|), 

A(b + f,bf) 



C 



A(B + F,bf) ' 



where the function A is defined in (10.9|. We then use the property (10.10) of A together with the 
estimate B + F < b + f from LemmajTl] to obtain that bf ^ 1 if and only if C ^ 1 . 



3.2 



we have B, F < 1. 



Finally, in Si?- interactions /, b < 1, and by parts (i) and (ii) o f Th eorem 
Again, using the explicit expressions for and </>, together with (4.3), we rewrite (|4Tj) as 

D(f) 



C = 



D(F) 



(4.9) 



where the function D is defined in (10.11 1 . Since D is strictly increasing and / > F by Lemma 



we obtain that C > 1 in SR interactions. This concludes the proof of part (iii) of Theorem 



3.2 



4.5. Proof of Theorem 3.2 part (iv). We argue as in Section 2.3 and observe that the map 



B i— > Q(B;b,f) is strictly increasing. The interaction Riemann problem is vacuum-free if and only 
if B = B(b, /) > 0, or equivalents G(B; b, f) = > £(0; 6, /), i.e. 



i>{b) + vM(b) > 



N(f) 



(4.10) 



Using the properties of the auxiliary functions ip, ip and M (see Section 10) we have 



for la and lb interactions, / < 1 such that (4.10) is satisfied 



LHS (4.10) > ip(b) > -v > 4){f) -v> RHS (4.10); 



• for Ic interactions b < 1 < / and (4.10) reduces to the condition that 

b^ + f~t > 1 (no vacuum) . 

The "vacuum curve" b^ + = 1 is sketched in the right diagram in Figure 2. This establishes 



part (iv) of Theorem 3.2 



This concludes the proof of Theorem |3.2| 

5. Proof of Theorem 13.31 



5.1. Group II: interactions with a contact. Consider the interactions in Group II listed in (3.2 ) 



Referring to the middle diagram in Figure 1 we use ( 2.6 )-( 2.8 ) to traverse the waves before and after 



interaction. We obtain three equations for the outgoing strengths B, C, F in terms of the incoming 
strengths / and c: 



c<K/) - C^(B)cb(F) 



(5.1) 

(5.2) 
(5.3) 



1>{f) = -^{B) + <ip{F)^CB4>{B) 
f = BF. 

Eliminating F and C, using the relation ( 10. 1[ ) and the definition (10.3), and rearranging, we obtain 
the following nonlinear equation for B = B(f, c) 

H(B;f,c)=0 (5.4) 



PAIRWISE WAVE INTERACTIONS IN IDEAL POLYTROPIC GASES 



15 



where 



f 



(5.5) 



We observe that since i\) is strictly increasing, the map B h-> T-l(B; f,c) has the same property. 
We proceed to analyze the strengths of the outgoing backward and forward waves B and F. This 



information is then used to verify the claims about the outgoing contact C in Theorem 3.3 We 



then establish part (iv) of Theorem 3.3 concerning vacuum formation in interactions with a contact. 
Finally, in Section 5.6 we include a partial analysis of the location of the transition curve {C = 1} 
in the (/, c)-plane, and show that it is the graph of a non-monotone function c n- f(c) satisfying 



limf(c) = limf(c) 

ell ctoo 



(5.6) 



5.2. Proof of Theorem 3.3 part (i). By the monotonicity of B i— >• H(B; /, c), the defining equa- 
tion (5.4) for B = B(f,c), and the definition (10.3) of the function N, we have 



B ^ 1 <=> 



0^K(l;/,c) = (l- v / H)^(/). 



Then, using the sign properties of ip, we obtain: f or S J interactions we have / < 1 and (3.7) foll ows, 
while for an RJ interactions we have / > 1, and (3.8) follows. This verifies part (i) of Theorem ° " 



3.3 



5.3. Proof of Theorem |3.3| part (ii). We observe that U(f; f, c) = V(/) + ip(f)- Thus, by ([5_3], 
the monotonicity of H(-;f,c), and the defining equation (5.4) for B = B(f,c), we have 



The last equivalence follows from the properties of the auxiliary functions ip and tp. This verifies 



part (ii) of Theorem 3.3 



5.4. Proof of Theorem 3.3 part (iii). As for Group I, in order to determine the type of the 



transmitted contact we will consider each interaction separately. The outgoing ratio C occurs in 



both (5.1) and (5.2), and we will make use of both relations. 



Case Ila. (S J) In this case /, c < 1, and from the analysis above we have that F < 1 < B. 
Thus / = BF > F, while a direct evaluation shows that B<fi(B) > 1. As is increasing we have 
> i>{f) > i>(F), and it follows that 



Mf) + 4>{B) > i>{F) > iP{F)^B<j){B) 



Comparing with (5.2) shows that C < 1 in this case. We further claim that C > c, which, according 



to (5.1 ), is the case if and only if 



cj>{f)=ct>{BF) > <j>{B)^(F). 



which is equivalent to E(BF) > E(B)E(F), where the function E is defined in (10.12). Wc apply 
part (a) of Lemma 10.2 and conclude that this last inequality is indeed satisfied since BF — f < 1. 
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Case lib. (S J) In this case / < 1 < c, and from the analysis above we have that F, B < 1. Thus 
/ = BF < F and the properties of the auxiliary functions tp an d show that 



A comparison with (5.2) shows that C > 1. We further claim that C < c, which, according to (5.1 1, 
is the case if and only if 

~$( f) = ~$ (BF) < %{B )${F) . 
Using the functions D, E (defined in ( |10.11 ) and (10.12)) we rewrite this as 

E{B)E(F) 



E(BF) < 



D(B) 



(5.7) 



Now, as B, F < 1, part (b) of Lemma [Kh2| shows that E(BF) < E(B)E(F). At the same time D 
is a strictly increasing function with D(l) = 1, such that D(B) < 1. It follows that (5.7) is satisfied 
and C < c. 

Case lie. (RJ) In this case c < 1 < /, whence, by the analysis above, B < 1 < F. This means 
that the backward and forward outgoing waves are both rarefactions, and it follows from Remark 
12.11 that C = c in this case. 



Case lid. (RJ) In this case c, / > 1, a nd the earlier a nalys is shows that 1 < F, B < f. Using the 
explicit expressions for <p and <p m (5.1), together with (5.3), we get 

C= J {S 1 C = ° (5.8) 
<t>(F)<f>(B) D(B) 

where D is defined in ( jlO.ll ). Since B > 1, D(l) = 1, and D is strictly increasing, it follows from 



(5.i 



that C < c in lid (RJ) interactions. 



Next we want to verify that C may be cither ^ 1 depending on the incoming strengths / and c. 



First, as B, F > 1 and BF = /, it follows from the defining relation %(B\ /, c) 
expressions for N , ip and ip, that 



(VE-i)/ c + i = v^B c 



and the explicit 



(5.9) 



The right-hand side is strictly increasing with respect to B, such that this relation provides B as a 
function of / and c: 



B = B(f, c) = unique B- value satisfying (5.9) for given /, c > 1. 



By substituting B = B(f,c) into (5.8) we obtain the strength C — C(f,c) of the outgoing contact. 
We now have: 

_^ < 

Lemma 5.1. Consider an RJ interaction with incoming strengths /, c > 1. Then, as f increases 
from 1 to oo while c > 1 is kept fixed, the strength C = C(f,c) of the outgoing contact (determined 
from (5.9) and (5.8)) decreases from c to 0. 

Proof. Recall that B i— > D(B) increases from 1 to oo as B increases from 1 to oo. Thus the claim 
follows from the expression (5.8) for C provided the map / B(f,c), for fixed c > 1, increases 
from 1 to oo as / increases from 1 to oo. Now, B(f,c) is the unique root of (5.9). For fixed c > 1 
the left-hand side of (5.9) increases from 1 to oo as / increases from 1 to oo. At the same time 
the right-hand side of (5.9) increases from 1 to oo as B increases from 1 to oo. This shows that 
/ i y B(f,c) increases from 1 to oo as / increases from 1 to oo. □ 



This establishes part (iii) of Theorem 3.3 
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5.5. Proof of Theorem 3.3 part (iv). We argue as in Section 2.3 and observe that the map 



B i y ~H(B; f, c) is strictly increasing. The interaction Riemann problem is vacuum-free if and only 
if B = B(f, c) > 0, or equivalent^ H(B; f,c) = > H(0; /, c), i.e. 

(5.10) 
Thus no 



u[l + y/iN(f)]-^(f)>0. 



Ila and lib: for an SJ interaction we have / < 1, such that LHS(5.10) > 0. 
vacuum occurs in S J interactions. 

lie and lid: for RJ interactions we have / > 1. Using the explicit expressions for the 
functions ip, ip, N, we get that no vacuum occurs if and only if 

(1-Vc)f ( <2. (5.11) 



For lid interactions c > 1, (5.11 1 is satisfied, and no vacuum occurs. On the other hand, for 



lie interactions we have c < 1, and no vacuum appears if and only if / < /*(c), where /*(c) 



is defined in (3.101 



This concludes the proof of Theorem |3.3| 

5.6. Location of the transition curve {C = 1} for lid interactions. We now want to determine 
more precisely the "transition" curve in the (/, c)-plane across which the outgoing contact changes 
type lid interactions. As above we denote the outgoing strength of the contact by C(f, c). By 
Lemma 5.1 we know that for each c > 1 there is a unique /-value f(c) such that C(f(c), c) = 1. We 
derive an expression for f(c) as follows. First, from ( |5.8[ ) with C = 1, we get that 

c = D(B). 

Let 5 := D^ 1 such that 

B = <5(c) = outgoing backward strength when no outgoing contact occurs. 



Substituting this value for B into (5.9 1 then yields the corresponding value of the incoming forward 
wave /, i.e. f(c): 

j; (6(c)) + tp (5(c)) 



f(c) 



5(c 



Kv^-i) 



where we have made use of the explicit expressions for ip. It turns out that this is a non-monotone 
function. Figure 7 shows the graph of the function 



T(B) := B c 



j,(B) + ^(B) 
u(^DjB)-l) 



when 7=|- (The plot shows the same qualitative features for other values of 7.) As c n- 5(c) is 
strictly increasing it follows that c 1— > f(c) is also non-monotone. In particular, the plot indicates 
that f(c) has a unique minimum. Finally, a direct evaluation verifies (5.6). 



Remark 5.2. We have not proved that f(c) has a unique minimum, nor determined its absolute 
minimum. In particular, we have not determined whether its absolute minimum is smaller or larger 
than f*(0) = 2 1 / c -. Figure 3 should be read with these provisions in mind. 



6. Proof of Theorem 3.4 setup and part (i) 



The analysis of interactions involving overtaking waves is more involved. From Theorems 3.2 
and |3.3| we see that most outcomes in Groups I and II interactions are independent of the adiabatic 
constant 7. Indeed, only the transition curve for the outgoing contact C in lid interactions, as well 
as the transition curves for vacuum in Ic and lie interactions, depend explicitly on the adiabatic 
constant. Also, these dependencies are "stable" in the sense that the transition curves are present 
and qualitatively similar for all values of 7 > 1. The situation for Group III interactions is markedly 
different. 



18 



CENG CHEN, ERIK E. ENDRES. AND HELGE KRISTIAN JENSSEN 



3.9 




100 200 300 400 500 600 700 SOO 



Figure 7. Graph of T(B) (7 = f ) 



We consider the pairwise interactions of overtaking backward waves of strengths x and y (see 
right diagram in Figure 1). It will turn out that the location and properties of the transition curve 
for the reflected (forward) wave, i.e. the locus {F(x,y) = 1}, depends sensitively on 7. At the same 
time, the transition curve of the outgoing backward wave, viz. {B(x,y) — 1}, as well as the vacuum 
transition curve, both depend on 7. However, in the latter cases the dependence is stable in the 
above sense. The situation is depicted in Figures 5 and 6. 

For reference we introduce the following three regions delimited by the lines x — 1 and y = 1 in 
the first quarter of the (x, y)-plane: 

• Ilia = {x, y > 1}, corresponding to S S'-interactions 

• Mb = {x > 1 > y > 0}, corresponding to S'i?-interactions 



• IIIc = {y > 1 > x > 0}, corresponding to .RS'-interactions. 

(The region {x, y < 1} corresponds to two backward rarefaction waves, which do not meet.) 

In this section we determine the type of transmitted (backward) wave in Group III interactions. 
Sections [7] and [8] provide the analysis of the reflected (forward) and contact waves, respectively, 
while Section [9] gives the conditions for vacuum formation in overtaking interactions. Together these 



results establish parts (i)-(iv) of Theorem 3.4 



6.1. Group III: overtaking interactions. Consider the interactions in Group III listed in (3.3 1 



Referring to the right diagram in Figure 1 we use (2.6)- (2.8) to traverse the waves before and after 



interaction. We obtain three equations for the outgoing strengths B, C, F in terms of the incoming 
strengths x and y: 



4>(x)cf>(y) = C<j>(B)<f>(F) 

^{x) + ^{y)^x%{x) = ^{B)-^{F)^JCB,AD) 
xy = BF . 



(6.1) 

(6.2) 
(6.3) 



Eliminating F and C, using ( 10.1 ) and ( 10.2 1, and rearranging, yield the following nonlinear equation 
for B = B(x,y) 

IC(B;x,y)=0, (6.4) 
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where 

K(B; x,y) := Tp(B) + % (M.) M (x) M (y) - ${x) - ^(y)M(x) . (6.5) 
We observe that since ip is strictly increasing, the map B H > IC(B] x, y) has the same property. 



6.2. Proof of Theorem 3.4 part (i). By the monotonicity of B H> 1C(B; x,y), the defining equa- 
tion (6.4) for B — B(x,y), and the definition (10.2) of the function M, we have 

B^l & 0^K(l;x,y) = M(x)K(x,y), 

where 

A(x,y) := i>(±)M(y) + - ^(y) . (6.6) 

We are not able in general to solve the equation K(x, y) = explicitly for x or y. Instead we 
determine the properties of the zero-level of K by estimating the partials d x K , d y K in Proposition 



6.1 We also determine the relative locations and intersections of {K(x, y) — 0} with the hyperbola 



xy = 1, see Proposition |6.2| This information will be used in Section [7] and depends on 7. 
Proposition 6.1. The partials of the function A(x,y) defined in \6.€fy satisfy 

d x K(x,y)<0, d y K(x,y)<0 for all x,y > . (6.7) 

Proof. As ip: "0 are strictly increasing we have 

d x K(x, y) = -^'(i)M(i,) - 4^'(i) < . (6.8) 

On the other hand 

9 B #(z, y) = -^'(^)M(y) + 4(±)M'(y) - i>'{y) , (6.9) 
such that d y K(x, y) < if and only if 

where the function to is analyzed in Section [10| By the properties of the auxiliary functions the last 
inequality is trivially satisfied whenever xy > 1. For xy < 1 we use instead that ip' > such that 
(6.10) follows provided 

xy^ (a;y) ^* ^ ( a:j/ / ^(2/) ' 

For xy < 1 the latter inequality is equivalent to 



where the auxiliary function I is analyzed in Section 10 In particular, £ satisfies £(q) > | for q > 1, 
while TO,(y) < I for all y > 0. This shows that d y K(x, y) < for all x, y > 0. □ 

It follows that {B = 1} = {A(x,y) = 0} is given by a graph y = fc(x), for a C 2 -smooth and 
strictly decreasing function fc(x). A calculation shows that: 

lim K(x,y) = 00, lim A(x,y) = —00 for x > 0, lim A(x,y) = 00 for y > 0, 

ylO vt°° X 10 

and 

lim K(x,y) = K^y) := -uM(y) - i()(y) 7= for y > 0. 

x-fco y a 

The function i^oo(y) is a strictly decreasing function and satisfies 

(0) = v 7= > , lim A'oo (y) = -00 , 



and has a unique root y = 2/(7) S (0, 1) (given in (9.5)). This concludes the proof of part (i) of 
Theorem [331 




Before considering the outgoing forward wave F we need to analyze the relative positions and 
intersections of the curves {xy = 1} and {B = 1} = {K(x,y) = 0}. These de pend on 7 and are 
given by the properties of the function a :— ip ~ i> which is analyzed in Section |l0. 3 



Proposition 6.2. Let {y = k(x)} = {B = 1} be as in part (i) of Theorem 3.4 Let y — 2/0(7) 
denote the root different from 1 of the function a(y, a) when 7 ^ |, and set yo = 1 when 7 = | . Let 
Xq = —. The relative positions and intersections of the curves {xy = 1} and {y = k(y)} are then 
given as follows: 

(a) For 1 < 7 < I (see Figure 8): yo < 1 < xq and 

• the curves intersect at (1, 1) tangentially and at (xo,yo) transver sally, and only at these 
points 

• \ < k(x) for x G (x ,oo) 

• i >k(x) forxe (0,l)U(l,z ) 

(b) For 7=| (see Figure 9): yo = x = 1 and 

• the curves intersect only at (1, 1) (tangentially) 

• \ < k(x) for x E (1, 00) 

• i > k(x) for x € (0, 1) 

(c) For 7 > § (see Figure 10): xq < 1 < yo and 

• the curves intersect at (1, 1) tangentially and at (xo,yo) transver sally, and only at these 
points 

• i < k(x) for x £ (x , 1) U (1, 00) 

• i > k(x) for x E (0,x ). 



Proof. From (|6.7|) we have that K(x, y) 



^ for y 5= 



k(x), while 
a(y,a) . 



(6.11) 



The statements about the location of the intersection points and the relative positions of {xy = 1} 
and {y = k(x )}, follow immediately from this and the properties of the function a(y, a) as detailed in 
Section 
that k'(l 



10.3 



q K ^ x ^ at the points of intersection shows 



(see Figures 8, 9, 10). Calculating fc'(x) = ~ 9 * K( - X > y 1 

^ — \ if and only if a' (yo) *s 0. The tangency and transversality 
claims follow from this and the properties of a. □ 
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rr^-r y(T) 



Figure 9. The sets {B = 1}, {xy — 1}, and their intersection for 7 = | (schematic). 




FIGURE 10. The sets {B = 1}, {xy — 1}, and their intersections for 7 > | (schematic) 



7. Proof of Theorem 3.4 part (ii) 
We now consider the reflected wave in overtaking interactions. To determine the type of the 



outgoing forward wave F we use (6.3) together with the defining equation (6.4 1 for B: 



(7.1) 
(7.2) 



F^l «■ xy^B H{x,y)^IC{B;x,y) = 0, 

where we have introduced the function 

H(x,y) := K(xy;x,y) = i/j(xy) - tp(x) - ip(y)M(x) . 
We note that 

H(l, y) = H(x, 1) = H(0, y) = Vx,y>0. 

To keep the leng ths of the proofs to a reasonable length we have collected some parts of the analysis 
10.5|l0.7 We consider separately the cases 7S(1,|), 7=f, and 7 > | . 



of H in Sections 



7.1. The reflected wave F in the case 1 < 7 < |. The following proposition details the proper- 
ties of the reflected wave when the adiabatic constant is between 1 and §. We recall that in := — . 

3 u yo ' 

where < yo < 1 is the unique root different from 1 of the function a = ip — tp defined in Section 

MM 
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F=l: 




B-l: 


yfi) h^is, 


,^("o,y ) 


; ^Z-——^ ~~~~~~ vacuum 



y.M 



Figure 11. The sets {F — 1}, {B — 1}, and their intersection for 7 < | (schematic) 



Proposition 7.1. Consider the interactions of two overtaking backward waves listed in (3.3). Let 
the left and right incoming waves have strengths x and y, respectively. For 1 < 7 < | (0 < a < j) 
the outgoing reflected wave F = F(x,y) is given as follows. 

(a) S S -interactions ('Ilia, x, y > 1) yield F > 1: the reflected wave is a rarefaction. 

(b) SB-interactions (TUb. < y < 1 < x) may yield either type of reflected wave. More 
precisely, in the region < y < 1 < x, the set {H(x,y) = 0} = {F = 1} coincides with a 
graph y = h{x), and the reflected wave F is a 

(bl) rarefaction (i.e. F > 1 ) if and only if y < h{x) 
(b2) shock (i.e. F < I) if and only if y > h{x). 

The location of the graph y = h(x) is given as follows. Let y — k{x) be the graph along which 
B = 1 (defined in Section \6.^ . Then, in the region < y < 1 < x: 

(b3) the three graphs y = h(x), y = k(x), and y — - all pass through (xo,yo), 

(b4) h(x) < k(x) < i for 1 < x < x 

(b5) h(x) > k(x) > \ for x > x a . 

(b6) y = h(x) has horizontal asymptote y*(^j) G (2/(7), 1) as x t +00, and 



lim h(x) 

2:^1 



4(l-a) 



< 1 



(7.3) 



(c) RS interactions (ITLc, § < x < \ < y) yield F < 1: the reflected wave is a shock. 
The situation is summarized in Figure 5 (left diagram) and in Figure 11. 



Proof. We consider each region in turn: 



(a) This part follows directly from (7.1) and (10.361 

(b) For this part we split the argument into several steps 

1 



We first observe that by definition (|10.2|) of the auxiliary function M, we have 
H( 



,y)=M(\)a{y,a) 



(7.4) 



Together with (6.3), ( |6.11[ ), and the properties of the map a(-,a) (see Figure 15), this 
shows: in the region Illb (0 < y < 1 < x), and for 1 < 7 < |, the two curves {xy = 1} 
and {B = 1} = {y = k(x)}, and the set {F = 1} = {H (x, y) = 0} (yet to be shown to 
be a curve), intersect at, and only at, the point (xo,yo)- 
Next we consider the equation H(x, y) = in the lower sub-region 

Illbl := {0 < y < xy < 1 < x} . 
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Using the explicit expressions for the auxiliary functions ip and M we may solve explic- 
itly for y as a function of x and obtain that: in Illbl the set {H(x, y) — 0} coincides 
with the graph y — h\(x) where 



hi(x) 



M(x) - $(x) - 1 
M{x) - scC 



\/x + a - *(x - 1) 



\/x + . 



X^y/jT 



(7.5) 



Note that the numerator and denominator are both positive by Lemma |10.1| A calcu- 
lation shows that h\(x) < ^ if and only if a(-, a) < 0, which, according to the analysis 
in Section 10.3 holds if and only if x < xq = ^ (since 7 < |). A direct evaluation 



shows that 



lim hi (x) 



4(l-a) 



yih) 



which is < 1 since a < \. This verifies (7.3). We conclude that {F = l}n Illbl is given 
as the graph {(x,hi(x)) \ 1 < x < xq}- Furthermore, it follows from ( |10.4 ) and (10.5) 
in Lemma TlO. II that 

{F ^ 1} n Illbl = {y§ hi{x), 1< x < x } . 

By the intersection properties verified in step 1 above, and the fact that 

lim k{x) = 1 , 

it follows that the graph {y = h\{x) \ 1 < x < xq} lies strictly below the curve segment 
{y = k(x) 1 1 < x < x Q }. 
3. The analysis of the equation H [x, y) = in the upper sub-region 

IIIb2 := {0 < y < 1 < xy} 

is more involved. In particular it is not possible to solve for y explicitly in terms of x 
or vice versa. Instead we will make use of the analysis in Section 10.6 First, (10.44) 
shows that the map y >— > H(x,y) has no root in (i, 1), whenever x < xq and 1 < 7 < |. 
We infer that {F = 1} n IIIb2 lies in the region x > xq. We next locate the zero-level 
of H more precisely. By ( |7.2[ ) we have H(x, y) — if and only if 

ip(xy) - tp(x) 



M{x) 



For (x,y) G {F = l}nlllb2 we thus get by ( |10.2| ) (see ([676]) for the definition of K) 



■ip (xy)- iji (x) 
M(x) 



V 1 (^v) 



< 0. 



cy/ ' w > M(x) 

That is, in sub-region IIIb2 the zero level of H (yet to be shown to be a graph y = h,2(x)) 
lies above y — k(x), and therefore (accor ding t o part (a) of Proposition 6.2) also above 
y = i. Finally, since 1 < 7 < f, ( 10.44| and |l0.45| show that {F(x,y) = 1} n {x > 
Xq, k{x) < y < 1} is a graph of a C^-function h,2(x). The graph y = h(x) in the 
statement of Proposition |7.1| is then the concatenation of y = h\(x) for x < xq with 
y = h%{x) for x > xq. 

4. We finally note that the inequalities between k(x) and =• were established in Part (a) 
of P roposition 6.2 This concludes the proof of part (b). 
(c) By (7.1) the reflected wave satisfies F < 1 if and only if H(x,y) < 0. When 7 < | we have 
a < j, and Lemma 10.5 shows that H [x, y) < in IIIc. 

□ 



Remark 7.2. While numerical plots indicate that the function y = h[x) is monotone decreasing, 
we have not been able to prove this. 
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7.2. The reflected wave F in the case 7 = |. This is a limiting case; for clarity of exposition 
we treat it separately. The following proposition details the properties of the reflected wave in the 
particular case when the adiabatic constant is that of a monatomic gas, 7 = S . 



Proposition 7.3. Consider the interactions of two overtaking backward waves listed in (3.3). Let 
the left and right incoming waves have strengths x and y, respectively. For 7 = | (a = 4 J the 
outgoing reflected wave F = F(x,y) is given as follows. 

(a) S S -interactions (Tlla, x, y > 1) yield F > 1; the reflected wave is a rarefaction. 

(b) S R-interactions (Tllb, < y < 1 < x) may yield either type of reflected wave. More 
precisely, in the region < y < 1 < x, the set {H(x,y) = 0} = {F — 1} coincides with a 
graph y = h{x), and the reflected wave F is a 

(bl) rarefaction (i.e. F > 1 ) if and only if y < h(x) 
(b2) shock (i.e. F < I) if and only if y > h{x). 

The location of the graph y = h(x) is given as follows. Let y = k{x) be the graph along which 

B = 1 (defined in Section \6.2\ . Then, in the region < y < 1 < x: 

(b3) the three graphs y = h(x), y = k{x), and y = - all pass through (1, 1), and 

(b4) h(x) > k(x) > i for x > 1. 

(c) RS interactions (lllc, < x < 1 < y) yield F < 1: the reflected wave is a shock. 
The situation is summarized in Figure 5 (right diagram) and in Figure 12. 

Proof. We consider each region in turn: 

(a) This part follows directly from (7.1) and ( 10.36[ ). 

(b) The proof of this part follows the proof for part (b) of Proposition 7.1 We split the argument 
into similar step s. 

1. From | |7.4| , (6.11), and the properties of the map ot{-\\) (see Figure 16) it follows 
that: in region Illb (0 < y < 1 < x) when 7 = |, the two curves {xy = 1} and 
{B = 1} = {y = k(x)}, and the set {F = 1} = {H(x, y) = 0} (yet to be shown to be a 
curve), do not intersect. 

2. Next we will show that the equation H(x, y) — has no solution in the lower sub-region 

Illbl := {0 < y < xy < 1 < x} 

whenever a > \. (W e will use this again in the proof of Proposition 7.4 ). As in the proof 
of Proposition I7T] we get that if (x, y) € Illbl satisfies H(x,y) = 0, then y 1 * = h\(x) 
where h\ is given by (|7.5[) (with a — j). In particular, since y < * in Illbl, we would 
have that h\{x) < -, which is equivalent to 

For a > j, or equivalently, 7 > |, this implies x < 1 (see Figures 16 and 17), which 
contradicts the assumption that (x,y) E Illbl. This shows that the set {F = 1} does 
not meet the sub-region Illbl when 7 > |. 

3. Next consider the equation H (x, y) — in the upper sub-region 

IIIb2 := {0 < y < 1< xy} . 
As in the proof of Proposition 



7.1 



10.6 



For a = I it is shown 



we use the results of Section 
there that the map y i-> H (x, y) has exactly one root in (^, 1) for each fixed x > 1 (recall 
.t = 1 in the present case). It follows that the set {F(x, y) — l}n{a; > 1, k(x) < y < 1} 
consists of a C 2 -smooth graph {y = h2{x)}. 

We finally note that the inequalities between A; (a;) and - were established in Part (b) 
of Proposition 6.2 This concludes the proof of part (b). 



(c) As in the proof of Proposition 7.1 this part follows from Lemma 10.5 



□ 




Figure 12. The sets {F — 1}, {B — 1}, and their intersection for 7 = | (schematic). 

7.3. The reflected wave F in the case 7 > |. This case is the most complicated one to analyze. 
In particular, now the set {F = 1} = {H = 0} either meets all three of Ilia, Illb, and IIIc, or just 
Ilia and IIIc. More precisely we shall see that {H(x,y) — 0} is a curve in the (x, y) plane which 
always meets Ilia and IIIc, and which meets Illb if and only if | < 7 < 2. See Figures 13 and 14. 

Proposition 7.4. Consider the interactions of two overtaking backward waves listed in t3Jm . Let 
the left and right incoming waves have strengths x and y, respectively. For 7 > | (a > j) the 
outgoing reflected wave F — F(x,y) is given as follows. 

(a) S S -interactions (Ilia,, x, y > 1) may yield either type of reflected wave. More precisely, 
in the region x, y > I, the set {F — 1} coincides with the graph of the strictly decreasing 
function 

2a 



Set 



j( x ) = Ti To- a(l + x) + \/a 2 (l + x) 2 + axil - a) 2 . (7.6) 

(1 — a) z x L J 

4a 2 _ 2a^L 4a 2 

Then the graph intersects the line x = 1 at (l,y). For ^<^<2(^<a<^) the graph 
intersects the line y = 1 at (x,l), while for ^ > 2 (a > ^) it has the horizontal asymptote 
y = > 1 as x j" 00 . The reflected wave is a: 
(al) rarefaction (F > 1) if and only if y > max(j(a;), 1) 
(a2) shock (F < \) if and only if 1 < y < j(x). 
(b) S R-interactions (Illb, < y < 1 < x) may or may not yield either type of reflected wave. 
First, for 7 > 2 (a > i) the set {F = 1} does not meet Illb and the reflected wave is 
necessarily a rarefaction (F > I). 

On the other hand, for |<7<2 (-g < a < ^) the set {F — 1} meets Illb along a graph 
y = h{x), defined for x > x, with h(x) = 1. The reflected wave F is a: 
(bl) rarefaction (F > I) if and only if y < min(/i(x), 1) 
(b2) shock (F < 1 ) if and only if x > x and h(x) < y < 1. 

The graph y — h{x) in the region < y < 1 < x lies above the graph y = k(x) along which 



B = 1 (defined in Section 6.2) 



(c) RS -interactions (lllc, < x < 1 < y) may yield either type of reflected wave. More 
precisely, in IIIc the set {F = 1} consists of a graph y — i{x), < x < I, with the 
properties: 
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(cl) a rarefaction results (F > 1) if and only if y < i(x) 
(c2) a shock results (F < 1 ) if and only if y > i{x) 
(c3) i(x) = yo (the unique zero of a(y, a) ) for < x < ^ 
(c4) i < i{x) <min{f,j/ } for i < x < 1. 
The situation is summarized in Figures 6, 13, and 14- 

Proof. We consider each region in turn: 

(a) Ilia: The set {F = 1} is the zero-level of the function H{x,y) whose behavior on Ilia is 
analyzed in Section 10.5 The expression in (7.6) is what results from solving (10.35) (with 
equality) for y in terms of x. A calculation shows that this is a decreasing function of 
x > 1, and that it has the intersection and asymptotic properties as described above. The 



conclusions follows from this together with (7.1) and (10.35 



(b) Illb: As shown in the proof of part (b) of Proposition 7.3 (step 2) the set {F = 1} does 
not meet subregion IIIbl:= {0 < y < xy < 1 < x} whenever a > 4. Also, according to the 
analysis in Section 10.6 {F = 1} meets subregion IIIb2:= {0 < y < 1 < x y < x} if and only 
The properties in (bl) and (b 2) follow from |l0.46[) and (|l0.47|. The fact 

kjxj 



if \ < a < k- 



4 

that y 
Proposition 7.1 (step 3) 



h(x) (denoted h2(x) in Section 10.6) lies above y 



is proved as in part (b) of 



(c) IIIc: It is convenient to consider separately the two sub-domains 

IIIcl := {0 < x < xy < 1 < y} and IIIc2 := {0 < x < 1 < xy < y} . 

In IIIcl we use the definition of H and the explicit expressions for the auxiliary functions ip 
and M to find that 



H(x, y) = ip(xy) - ip(x) - ip(y)M(x) = x c a(y, a) , 

where a := ip — tp has a unique zero yo > 1 (since 7 > |, see analysis in Section 
Figure 17). The properties of a then shows that the restriction of H to IIIcl satisfies 



10.3 



and 



H(x,y) 



< for < x < i and y > yo 



= along (x, y ), < x < 



(7. 



> for < x < i and 1 < y < y Q . 



In sub-region IIIc2 we consider instead how H(x,y) varies as (x,y) moves along hyperbolas 
{xy = C} (C — constant) in the direction of increasing j/-values. For this we make use of 
the properties of the function rj{y,a), which is analyzed in Section 10.7 In the rest of this 
proof we assume C > 1. Let's define the directional derivative 



f V) : = V) ■ ^(x,y)H(x, y) 



%x c rj(y,a) , 



where rj(y,a) is defined in (10.49). We first observe from ( 10.48 )-( 10.49 1 that the leading 
term in 0(— ,y) for y S> 1 is proportional to — y 3 ~^. Since C < | ^ ms snows that j/) | IIIc2 
tends to —00 as y f 00. 

We next consider the sign of £)(x, y) as we "start out" along {xy = C} from (1, C), in the 
direction of increasing y. As detailed in Section 10.7 the sign of 5(a;, y) coincides with that 
of 



1 - Va 



\Jy + a. 



Thus, if the constant C satisfies 
\/a(l + a) 



a<0 ^> C >y :-- 



2a^/a 

(i-VS) s 
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Figure 13. The sets {F = 1}, {B = 1}, and their intersection for | < 7 < 2 (schematic). 



then, since H(1,C) = 0, H(x,y) < along {xy — C}. On the other hand, if 1 < C < y, 
then H increases along {xy — C} as y increases, until y = y, after which it decreases to 
—00. We see from this that H has: 

— no zero along {xy — C} when C > y, 

— a unique zero along {xy = 0} when 1 < C < y. 

It follows that the zeros of H(x,y) in IIIc2 lie along a curve y = i(x), where i satisfies 
i(^) = y and = y. It remains to argue that the graph y — i[x) lies below the line 
y = j/o, and we do this by showing that H(x,yo) < for ^ < x < 1. Indeed, by using the 

property ip(yo) = i>(yo) together with the explicit expressions for the auxiliary functions tp 
and M, we have that 

d x H(x, y ) = -d x (a(xy 0) a)) for — < x < 1. 

yo 

Integrating from i to x we obtain 

H(x,y ) = -a(xy ,a) < 0, 
where the latter inequality follows from the properties of a(-,a) when a > 4. 

□ 



Proof of Theorem 3.4 part (iii) 



To analyze the outgoing contact discontinuity in Group III interactions it is advantageous to work 
in (t, u, S')-space where we track two quantities, r and S 1 , that change across contacts. The wave 
curves in these variables were recorded in Section [2j 

8.1. Equations for outgoing waves in (r, U, S')-variables. We refer to Group III interactions in 
Figure 1 and denote the specific volumes ratios T r i g ht/Tieft across the incoming waves by x (leftmost) 
and y (rightmost). We use capital letters L, C, I to denote the specific volume ratios across the 
outgoing backward, contact, and forward waves, respectively. 

We use the expressions for the wave curves in (t, u, S*)-space to traverse the waves before and 



after interaction. From (2.10)-(2.13) we obtain the following equations for the outgoing strengths 



28 



CENG 



CHEN, ERIK E. 



ENDRES, AND HELGE KRISTIAN JENSSEN 




Figure 14. The sets {F = 1}, {B = 1}, and their intersection for 7 > 2 (schematic). 



L, I, C: 



rj(x) + f)(y) = Tf(L)- 
xy = CLI . 



7 log C +rf(I), 



(8.1) 
(8.2) 
(8.3) 



As for Groups I and II we consider Ilia, Illb, and IIIc interactions separately. We find it necessary 
to make a further breakdown and consider each combination of outgoing extreme waves within each 
of Ili a, II Ib, and IIIc. For the most part the type of the outgoing contact follows readily from 
.1 )-(8.3 1 and the properties of the auxiliary functions. The only exception is the case S S —> SJR 



which requires additional arguments. 



8.2. Outgoing contact in Ilia-interactions. These are S'S'-interactions, which corresponds to 
a < x, y < 1, or, in terms of incoming pressure ratios x, y > 1. It follows from the analysis in Section 
6.2| that, independently of the value of 7, the outgoing backward wave is a shock, i.e. a < L < 1. 
There are therefore only two possibilities for the extreme outgoing waves in this case. (The analysis 
in Section [7] shows that both can occur when 7 > |.) We treat them separately and show that the 
outgoing contact satisfies C < 1 in both cases. 



Case 1: SJS. In this case a<L<l<I<\. By (8.2 1 and (10.25) we have that 



T{x)T{y) - T{L)T{I)& 
Assume for contradiction that C > 1; then 



(8.4) 



4) together with (10.26) give 

r(L) < mrimi- 1 ) ■ 



Also, if C > 1 then (8.3) gives xy > LI > a. Since a < x, y, xy < 1, Lemma 10.3 gives r(x)T(y) < 
T(xy), such that 

r(L) <r(^)r(/- 1 ). 



10.3 



Agai n, since C > 1, (8.3 1 gives xyl 1 = CL > L > a, and we have a < xy, I l , xyl 1 < 1. Lemma 
applies and gives r(iy)r(/ _1 ) < T^yl^ 1 ), such that 

r(L) <T(xyr l ). 



As T is strictly decreasing on (a, -), we conclude that L > 7^. However, by (8.3), this implies that 
C < 1 and we reach a contradiction. Thus C < 1. 
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Case 2: SJR. In this case a < L < 1 and I < 1. As x < 1 and £ (I) > 0, ( |8.l[ ) and ( |10.25[ ) give 
that 

(8.5) 



Z(x)+Z (y)Vr{£) < £ (5) + £ (y)v^f(I) < £(L). 



Also, in this case (10.25) shows that (8.2 1 reduces to 



(8.6) 



T(i)r(y) = T(L)C . 
Assuming, again for contradiction, that C > 1 we thus obtain 

T(x)T(y) > T(L) . 

We proceed to show that this leads to a contradiction with (8.5). For this we set f2(s) := v/r(s) 
(s G (a, 1)) and define the function A := £ (ri _1 ) which is analyzed in Section 10.4 Now let 

zx := 0(5) > 1, z 2 := ft(y) > 1, z 3 := > 1 , 



where we have used that T and hence are decreasing. Hence (8.6) reduces to 

z l z 2 > z 3 j 



(8.7) 



and Lemma f!0.4| gives 



A(2i) + ziA(z 2 ) > A(ziz 2 ). 



Thus, by (8.7) and the fact that A(z) is an increasing function, we obtain A(zi) + ziA(z 2 ) > A(zs), 
or equivalently: 

t(x) + VWH(y) > • (8.8) 



This contradicts (8.5), and we conclude that C < 1. 



This establishes the first statement in part (iii) of Theorem 3.4 



8.3. Outgoing contact in Illb-interactions. These are S'i?-interactions for which a < x < 1 < y. 
There are now four possible combinations of outgoing forward and backward waves. (The analysis 
in Section [7] shows that they can all occur when 7 > |). We demonstrate that the outgoing contact 
discontinuity C always satisfies C > 1 by considering each case separately. 



Case 1: RJS. In this case L > 1, 1 < I < - and (10.25) shows that (8.2) reduces to 



r(5) = r(/)c 7 



Since a < x < 1 < I, (10.281 shows that T(x) > T(I), such that C > 1. 



Case 2: RJR. In this case L > 1, I <1 and (10.25) shows that (8.2) reduces to 

r(i) = C 7 . 



Since a < x < 1, (10.28) shows that T(i) > 1, such that C > 1. 



Case 3: SJR. In this case a < £ < 1, I <1 and (10.25) shows that (8.2) reduces to 



F(x) = T(L)C 7 



We now argue by contradiction: if C < 1 then T(x) < T(L), and ( 10.28 ) shows that L < x. However, 
combining this with C < 1 < y and (8.3) gives 



L<Ly<xy = CLI < LI, 



contradicting I < 1. Hence C > 1. 



Case 4: 5 JS\ In this case a<L<l</<±and (8.2), (10.25), and (10.28) give 



T(x) = T(L)T(I)C < r(L)C 7 . 



(8.9) 



1) and the fact that £ (y), £(7) < 0, we obtain that £ (i) > £ (£). As £ is strictly 



^From 

decreasing, we infer that x < L. Thus T(i) > T(L), which combined with (|8.9[) gives C > 1 
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8.4. Outgoing contact in Illc-interactions. These are i?S'-interactions for which a < y < 1 < x. 
As for S'.R-interactions there are four possible combinations of extreme outgoing waves. (The analysis 
in Section [7] shows that they can all occur when 7 > |). We claim that the outgoing contact 
discontinuity C always satisfies C > 1. Considering the same cases as for S'i?-inter actions it turns 
out that the arguments for Cases 1, 2, and 3 for .RS'-interactions are identical to those for SR- 
interactions, upon interchanging x and y. We therefore only need to c onsi de r Case 4, S JS, wh ere 
the outgoing strengths satisfy a<L<l</<i. As above we use 
obtain 

r(y) = r(L)r(/)c 7 < r(L)c 7 . 

Since x > 1 we have £ (x) < and y'x 1-7 exp(rf(x)) < 1. At the same time | (J) < 0, and we get 
from 



5.2|, (10.251, and (10.28) to 



(8.10) 



.1 ) that £ (y) > £ (L). As £ is strictly decreasing we have y < L and (8.10) yields C > 1. 



With this we have established the second statement in part (iii) of Theorem 3.4 



9. Proof of Theorem 3.4 part (iv) 

For this part of the proof we use the pressure ratios x and y of the incoming waves. We argue as in 
Section 2.3 and observe that the map B h-> IC(B; x, y) is strictly increasing (see (6.5 ) for the definition 
of the function K,). The interaction Riemann problem is vacuum-free if and only if B = B(x, y) > 0, 
or equivalently K,(B;x,y) = > /C(0;x, y), i.e. 



tC{0;x,y) = — v — vM{x)M{y) - tp(x) - rp(y)M(x) < 0. 



(9.1) 



Rearranging the last expression we have that the overtaking-wave interaction produces no vacuum 
if and only if the incoming parameters x and y satisfy 



vM(y)-ip(y) < 



v + tp(x) 
M(x) 



(9.2) 



Sinc e the function ip is strictly increasing with ^(0) = — v and "0(1) = 0, it follows that the inequality 
(9.2) is satisfied whenever y > 1. Consequently, a vacuum is never generated in S S (Ilia) and RS 
(IIIc) interactions. 

On the other hand, depending on the incoming pressure ratios x a nd y , a vacuum may or may 
not emerge from an SR (Mb) interaction. In this case y < 1 < x and (9.2) takes the explicit form 

v\Jx + a + k(x — 1) 



v(x) 



Vx + 1 



> K1-2/) 



(9.3) 



Lemma 9.1. The function v(x) is strictly decreasing for x > 1, for all values of ^7 > 1. 



Proof. Differentiating and collecting terms with coefficients v and k, respectively, we obtain that 
v'(x) < (for x > 1) if and only if 

I x + a 1 + 2ax + x 2 



l + a < 1 + (1 + 2a)x ' 
This relation holds since the right- and left-hand sides are separated by a linear function: 

1 + 2ax + x 2 



x + a x — 1 
1 + a < 2(1 + a) 



1 < 



1 + (1 + 2a)x 



for x > 1 . 



□ 



Thus, a vacuum appears in a Mb interaction if and only if the incoming strengths y < 1 < x 
satisfy 



0<y<V(x) := [i (i-^i) 



where C 



7-1 
27 ' 



(9.4) 
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and v(x) is given by (9.3 1. A calculation shows that the strictly increasing function V(x) satisfies 
V(l) = and 

lim V(x) = [i (1 - VC)] * =: y( 7 ) < 1 . (9.5) 

xfoo 

See Figures 5, 6, 8, 9-14 for schematic plots of the vacuum transition-curve y = V{x). As indicated 
in these figures the horizontal asymptote y("f) of y = V{x) coincides with the horizontal asymptote 
of the transition curve {B = 1} in SR (Illb) interactions. 

This completes the proof of Theorem |3.4| 



10. Definitions and properties of auxiliary functions 

10.1. Auxiliary functions for wave curves. The functions 4>, <f>, ip, ip were defined in Section 
[2J A calculation shows that they are all C 2 functions with Lipschitz continuous 2nd derivatives. 
Furthermore: 

• 4>{q) is strictly decreasing, tends to +00 as q \, 0, tends to a as q — > 00, and </>(l) = 1. 

• (j)(q) is strictly decreasing, tends to i as q | 0, tends to as q — > 00, and </>(l) = 1. 

• ^p{q) is increasing, tends to — v (with infinite slope) as q 4- 0, tends to +00 as q — > 00, and 
5>(1) = 0. 

• ?p(q) is increasing, tends to — \J^f L (with finite slope) as q \. 0, tends to +00 as q — > 00, 
and ^(1) = 0. 

For reference we record the relation 

=-?(<?) • (10-1) 



10.2. The functions M, N, m, n, £, A, D and E, and their properties. In this section we 
define a number of auxiliary functions and list some useful properties. To verify these requires mostly 
routine calculations which are not included. We recall that the parameter a € (0, 1) is defined in 



2.9 Define the functions M and N by 



and 



M{q) 



N{q) 



HI) 



VKg) 

>0 



q+aq A 
q+a 



q+aq 2 
q+a 



0<q<l 
q> 1 

< <? < 1 
?>1 



(10.2) 



(10.3) 



Both M(g) and iV(g) take the value 1 at g = 1, are increasing and convex down, tend to zero as 
q X (with infinite slope), and tend to +00 as g — >• 00. 



Lemma 10.1. For a; > 1 the function M{x) defined in \10.2 ) satisfies the following inequalities ft 
all values of 7 > 1 : 



or 



M{x) 
M(x) 



'x -+ 


ax 2 


X 


+ a 


!x + 


ax 2 



x + a 



> x K 



> 1 



k(x — 1) 
v\Jx + a 



(10.4) 
(10.5) 
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Proof. By squaring both sides in ( 10.4 1 and rearranging we obtain the equivalent condition that the 
function 

f{x) :=x + ax 2 ~ x 2C+1 - ax 2C 

satisfies f(x) > for x > 1. A calculation shows that f"'(x) > for x > 1, while f(l) = f'(l) = 
/"(l) = 0. This proves (10.4). By squaring, rearranging, and canceling a factor of (x — 1) we obtain 
that ( 10.5 ) holds if and only if 



(l-C)* + (l + 0>2^±f 
Squaring again and simplifying gives that this holds if and only if (x — l) 2 > 0. □ 
Define the function m by 

„ x f =C < q < 1, 

M{q) 1 _ a(l+2aq+q 2 ) -, 
^ _ 2(q+a)(l+aq) 1 > 1 • 

A calculation shows that m is non-decreasing and with range [£, ^). Hence the function 

n{q):= ^w =m ^ {q>0) (10J) 

is non-increasing and with range [£, J). Define the function 



_£2l 



<<?< 1, 



«:=^= (10.8) 

1p{q) q(q+2a+l) , 

I 2( g -l)(g+a) q > 1 ' 

A calculation shows that ^ is decreasing, has a vertical asymptote at q — 1, and satisfies i(q) > | 
for all g > 1. Define the function A by 

A(q,0 := 1 ^ + " 2 f (g,f>0). (10.9) 

A calculation shows that 

d q A(q,0^0 Vq>0 e> e^l. (10.10) 

Define the function I? by 

D(g) := ( i^?) (g>0). (10.11) 

A calculation shows that D is strictly increasing, D(l) = 1, and that D(q) f oo as q f oo. 
Lemma 10.2. Define the function E by 

£( x ):=^^ x>0 . (10.12) 
x + a 

Then 

(a) /or < x < 1 < y: E(xy) > E{x)E(y) ^ xy < 1, 

(b) for < x, y < I: E(xy) < E(x)E(y). 

Proof. A calculation shows that, since < a < 1, E(xy) > E(x)E(y) if and only if (1 — xy)(l — 
x)(l — y) < 0. Parts (a) and (b) follow directly from this. □ 
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such that 



10.3. The function a and its properties. The function a plays a key role in several parts of the 
arguments. We define 

a{q) = a(q, a) := $(q) - ^(q) q > 0. (10.13) 

We need to locate the roots of a, and these depend sensitively on the value of a — -^j. As a first 
step we introduce 

a(q) = a(q,a):='^B-v( q t-l), q>0, (10.14) 
V q + a 

a(?) = | _ &{q) (10-15) 

where the constants k, v, ( are defined in terms of a in Section [2j To analyze a we introduce the 
new variable 

*(«) = *(9,a):=4=. q>0, (10.16) 
VI + a 

which satisfies 

,. q + 2a+l „ q + 4a + 3 /in ^\ 

z (l) = h nr, z («=-77 nr- (10-17) 

2(q + a)3 4(g + a)2 

For fixed a € (0, 1) the function z(q) is strictly increasing with a strictly increasing inverse q(z) 
defined for z € (— a -1 / 2 , oo), and with q(0) = 1. The first two derivatives of the latter are given by 



2 



q > (z) = J_ = 2 (g + fl ) 5 ^(,) = - " //(9) = % + 4a + 3)(9 + a) 

91 J z'(q) q + 2<z + l' 91 j (*'(q)) 3 (g + 2a + 1)3 ' [WAti) 

where q = q{z). Now, in terms of the variable z(q), we have 

a(q) = f3(z(q)) , 

where 

P(z) := nz-v[q{zf (10.19) 

Thus, in order to locate the zeros of a we may as well determine the zeros of (3(z) for z 6 (— a -1 / 2 , oo), 
and then translate back to q-locations. (This turns out to be easier than to determine directly the 
roots of a.) We do so by considering the derivatives of f3(z). Differentiating j3{z) and using (10.18) 
we have 

VTTa _ q^ 1 {q + a)i 
2~ q + 2a + 1 



P'(z) = 2K 
and 



(10.20) 



„„, N 2q c ~ 2 (<7 + a) 2 
75(q + 2a+ l) 3 

where 

_ 2a(2a + l) 
q = q(z) ana q :=■ 



1 -a 

This shows that the sign of (3"(z) is the same as that of the function 

Q{z) := (q(z) - l)(q(z) - q) , (10.21) 

such that Q(z) has the zeros z\ := and Z2 := z(q). We note that 

> ->i > I >5 nn , 

^2 = 2i q=l O a=- ^ 7=-. (10.22) 

To locate the zeros of j3(z) we consider these regimes separately. The approach is the same in each 
case: the signs of f3"(z), together with end-point values of f3" , (3' and j3, determine the number and 
locations of the roots of j3. We therefore only detail the argument in the representative case when 
1<7<§- 
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FIGURE 15. q h> a(q, a) when 0<a<j(l<7<|, schematic) 




Figure 16. q i-» a(q, \) (7 = |, schematic) 



In t his case the function Q(z) has two distinct roots z\ = and z 2 = z(q) G (—a 1 / 2 ,0). ^From 
( |l0.2l[ ) we have that /3"(z) > when z e (-a -1 / 2 ,^) U (0,oo), and /3"(z) < when z £ (z 2 ,0) . 
Hence f3'(z) is increasing on (— a _1 / 2 ,z 2 ) U (0, 00) and decreasing in (z2,0). Next, from (10.16 1, 
(10.191, and £ < ~ we obtain 



lim 



-00, 



/3'(0) = 0, lim p{z) = Kv^Ta > 0. 



(10.23) 

, Z2) and one at z = 0. Also, /3'(z) < on 



It follows that (3'{z) > has two roots: one at z 3 G (— < 
(-a _ 5,z 3 ), and (j'j z) > on (z 3 ,0) U (0,oo). 

Next, by jl0.16| ) and ( |l0.19| we have ^(-a^ 1 / 2 ) = 1/ - ko" 1 / 2 > and /?(0) = 0. It follows 
that j3 itself has exactly two roots: a leftmost root 24 € (— a -1 / 2 , 23) and the other one at z = 0. 
Translating back to q- variables we conclude that a(q) has exactly two zeros: (7(24) =: yo € (0, 1) and 
q=l. Finally, it follows from this and (10.15) that the same holds for a(q) itself. See Figure 15. 

Similar arguments show that 

• when 7=| the function a(q) has a single root at q — 1 (Figure 16), and 

• when 7 > I the function a(q) has exactly two roots: one at q = 1 and one at yo £ (l,oo) 
(Figure 17). 



PAIRWISE WAVE INTERACTIONS IN IDEAL POLYTROPIC GASES 



35 



a(q) 



Figure 17. q h-» a(g, a) when a > | (7 > |, schematic) 
10.4. The functions V and A, and their properties. We define the function 



r( s ) = r( s ,<z) 



sT(l - as) 



a < s < 



1 







such that the functions rf and 77 in (2.11 1-(2.13| may be expressed as follows: 
( 1>1 ) 

m = { \ -- 

[ io g r(/) a<i<i J [ io g r(t 

The function T satisfies the relation 




r( s )r(i) = i, 



1 

a < s < - 
a 



A calculation shows that 



r'( s ) 



-07s ■ 



(s - a) 2 ' 

such that r is a positive and strictly decreasing function on (a, and satisfies 

r(s) ^ 1 for s $ 1. 
Lemma 10.3. 27ie function T has the following property: for s, t, st 6 (a, A) 

r(s)r(t) < r(st) tf and only if (1 - st)(l - s)(l - t) > 0. 
Proof. By using the expression for T we have that r(s)r(i) < r(si) if and only if 



1 — as 



1 - at 

t — a 



< 



1 — ast 

st — a 



(10.24) 

(10.25) 

(10.26) 
(10.27) 

(10.28) 
(10.29) 



Using the assumption that s, t, and st are all larger than a, a calculation shows that this 
and only if (1 - st)(l - s)(l - f) > 0. 

We introduce the auxiliary function 

A(z) := *€(l,oo), 

where £ is defined in Section 



holds if 
□ 



(10.30) 



2.2 



and 

n(s) := y/T(sj, sg(o,1). 
Lemma 10.4. T/ie function A : (1, 00) — > (0, 00) is increasing and has the property that 
A(zi) + ZiA(#2) — -M-Zi^) /or a// 21, z 2 > 1 . 



(10.31) 
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Proof. First, since A(1)=0, (10.311 follows from 



A'(z) > A'(ziz), z e (l,oo) 



(10.32) 



by integration from z = 1 to z = Z2- In turn, ( 10.32 1 follows if we show that A is a concave: 

A"(z)<0, z e (l,oo). 



To establish (10.33) we first calculate A'(z) 



A'(z) = e ' (s) = Vl^i 8 + 1 - 2a ) 

Z> SV(a) 07 ' si- x (s-l) 2 



where s = 1 (z) 6 (a, 1) 



Since O'(s) < we see that (10.33) holds if and only if 

d 
ds 



Vl - as(s + 1 - 2a) 



> for s e (a, 1). 



s i-i(s- 1)2 

A calculation shows that this is the case if and only if the polynomial 

,2„3 /in„2 r„ 1 1 \„2 1 /o„3 in„2 , i ,4 „ r\ „ (c„1 



satisfies 



V(a, s) := 2aV - (10cT - 5a + l)s 2 + (8a - 10a^ + 14a - 6)s - (6a 2 - 5a + 1) 



■p(a, s) < whenever < a < 1 and a < s < 1 . 



(10.33) 



(10.34) 



To verify ( 10.34) we fix s g (0, 1) and consider the map a i-> V(a, s) for < a < 1. Since 
d 3 a T(a, s) = 48s > 0, d 2 a T{l, s) = 4(s - l) 2 (s - 3) < 0, 

while 

d a V{s, s) = (s- l) 2 (4s 2 + 12s + 5) > and T{s, s) = (s - l) 3 (2s 2 + 4s + 1) < , 



it follows that (10.34) indeed holds 



□ 



10.5. The function H(x, y) in region Ilia. Consider the function H defined in (7.2). For (x,y) £ 
Ilia = {a; > 1, y > 1} we use the explicit expressions for -0 and M to get 



H(x,y) = k 



(xy — 1) (x — 1) (y — 1) /x + ax 2 
^/xy + a ^Jx + a \/y + a V x + a 



(x, y > 1). 



Thus 



^xy + a < y/x + a ' ^/y + a V x + a 

(xy — l)\fx + a\Jy + a ^ (x - l)y^ + a-^/xy + a + (y - l)y/xy + a\J x + ax 2 . 

For x, y > 1 this holds if and only if the same inequality with squared RHS and LHS holds. After 
squaring each side, collecting positive and negative terms, canceling the common factor y / x(x — 
l)(y — 1), and rearranging, we obtain: 

H(x,y)^0 4=> 2(xy + a)y/y + a\A + ax ^ Jx [(1 - a)xy 2 + 2a(xy + a) + (3a + l)y] . 

Again each side is positive; squaring, collecting terms, canceling the common factor (1 — xy) 2 , and 
rearranging, finally yield 

7J(x,y)^0 ^ xy [(1 - a) 2 y - 4a 2 ] ^ 4a 2 (y + a) when x, y > 1. (10.35) 

In particular, if < a < \ then (1 — a) 2 y — 4a 2 > 4a 2 (y + a) whenever y > 1. Thus 

H(x,y) > whenever x, y > 1 and < a < j. (10.36) 
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10.6. The function H(x,y) in region IIIb2. Consider the function H defined in (7.2 1. We want 
to determine the its zeros in the subregion IIIb2 := {0 < y < 1 < xy}. ^From (7.2) and the explicit 
expressions for if) and M, we have 



H(x, y) = - - v(y< -1)J X + 



x + a 



for — < y < 1. 
x 



(10.37) 



yjxy + a y/x + a 

Fixing x > 1 we first note that 

H(x, i) = M(x)a(A,a) while ff(i,l) = 0, Vx>0. 

We study y i— > H(x,y) as y increases from ^ to 1 by analyzing y i— > d y H(x,y). As we shall see, its 
behavior depends on whether < a < |, | <a< |, or a > |. Introducing the functions 9(z, a) 
and Q(z, x, a) by 

(z + 2a + l);? 1 ^ 



Q(z,x,a) 



(z + a) 5 

V 7 ! + a x c 
2 M(x) ' 



for z > 1, 



(z, a) — 1 for x, z > 1, 



we have 



dyff (x, y) = vCM(x)y c - 1 Q{xy, x, a) for (x, y) € IIIb2. 
We next observe the following points: 



(10.38) 
(10.39) 
(10.40) 



(Al) By (10.40) the sign of d y H(x,y) is the same as that of Q(z, x, a)\ z 
(A2) According to Lemma 10. 1| 



Q(l I x ! a)= J ^-l<0, 

such that, in IIIb2, y \-> H(x,y) "starts out" decreasing at y = ^. 
(A3) A calculation shows that Q(x, x, a) > for x > 1 if and only if 



(10.41) 



(1 — 3a) (x — x) > where x 



4a 2 



l-3a 



(10.42) 



(A4) We have x < for a > |, while 

0<x<l ^ < a < j , x>l ^ ae(|,|). 

For a = \ , x is undefined and Q(x, x, |) < 0. 
(A5) By ( 10.39 1, the sign of d z Q(z, x, a) coincides with that of d z 8(z, a). A calculation shows that 

2a(2a + 1) 



the latter is given by 



d z 9(z,a) = 



(l-a)(z-l)(z~z) 
2(l + a)zC(z + a)i 



where 



1 - a 



(10.43) 



We have 



z^l 



1 



We can now analyze the zeros of y h- > H(x,y) for - < y < 1: 
• < a < |. By (A2)-(A5) we have 

Q(l,x, a) < and Q(x,x,a) > 0, while d z Q(z 1 x 1 a) > for 1 < z < x. 



It follows from this and (10.40) that y i-> d v H(x 1 y) changes sign once from negative to 
positive as y increases from ^ to 1. From (10.371 and the properties of the map a, we 
have H(x, ~) ^ according to i ^ x . (Here Xq = — , yo being the unique root of a(-,a) 
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different from 1 when a < j, and y n = 1 when a=\-) We conclude from (10.371 that when 
1 < 7 < | then the map y h- > fi(:r, y) has 



no root in (A, 1) when a: < xo, 
exactly one root y = /i2(x) € (-, 1) when x > xq. 

As H is a C 2 -map it follows that /i2(x) is a C 2 -function. 



(10.44) 
(10.45) 



• a > \. By (10.371 and the properties of a(y,a) we have H(x, A) > = fi(x, 1) for a; > 1. 
As above we have Q(l, x, a) < 0, but now z > 1 (by (A5)) and 

d z Q(z, x, a) < for z e (1, z), d z Q(z, x, a) > for z e (z, oo), 



and Q(z,x, a) — » oo as z t oo (by (10.391, (10.381, and the fact that ( < |). Thus, the map 
z i— > Q{z,x,a) has a unique zero in (z, oo) whenever a > j. By ( 10.40 1 this root corresponds 
to a zero of y H> d y H(x,y) in the interval (i, 1) if and only if Q(x,x,a) > 0. By (A4) and 
(10.42) this is the case if and only if a < i and a; > i. We conclude that, for each fixed 
x > 1 the map y i-> if (a;, y) has: 

• a unique root y = h,2(x) € (^, 1) if and only if x > 5, when | < a < |, (10.46) 

• no root in (^, 1) when a > ^. (10.47) 

Again, in the former case hi is C 2 -smooth. 



10.7. The function H(x,y) in region IIIc. Consider the function H defined in (7.2 1. In IIIc:= 

{0 < x < 1 < y} we analyze how H(x,y) varies along hyperbolas xy = cords. Again, by using the 
explicit expressions for ip and M we have 



{-x, y) ■ V^,y)H(x, y) = %x c r](y, a) for < x < 1 and y > 1, 
where the function 77 is given by 



»7(y»a) : = 



(1 - a)y 2 + (5a + l)y + 2c 



\/l + a (l + a)(j/ + a)i 

A calculation shows that r](y, a) factors as follows: 



for y > 1. 



(10.48) 
(10.49) 



v(y,a) 



(1 - a) (7^+^- Vl + a)' 
(l + a)(y + a)i 



Vy + « 



Vo(l + a) 
1 - 



+ « + 



Va(l + a) 
1 + V« 



It follows from this and ( 10.48 ) that the directional derivative (— x, y) ■ Vt x , y -\H(x, y) in IIIc has the 
same sign as 



y/a(l + a) 



which is strictly decreasing in y. In particular, since 
^a(l + a) 



1 - 



- y/y + a 



< if and only if a < 



4 ■ 



y=i 



and since H(l,y) — H(x, 1) = 0, we obtain the following: 

Lemma 10.5. For a < \ the function H(x,y) is strictly negative in IIIc = {0 < x < 1 < y}. 
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